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We perform a Chapman-Enskog expansion of the Boltzmann equation keeping up 
to quadratic contributions. We obtain a generalized nonlinear Kubo formula, and 
a set of integral equations which resum ladder and extended ladder diagrams. We 
show that these two equations have exactly the same structure, and thus provide a 
diagrammatic interpretation of the Chapman-Enskog expansion of the Boltzmann 
equation, up to quadratic order. 



Fluctuations occur in a system perturbed slightly away from equilibrium. 
The responses to these fluctuations are described by transport coefficients u. 
The investigation of transport coefficients in high temperature gauge theories 
is important in cosmological applications such as electroweak baryogenesis 
and in the context of heavy ion collisions . There are two basic methods to 

Sulate transport coefficients: transport theory and linear response theory 
I. To date, most calculations of transport coefficients have been done to 
the order of linear response. In many physical situations however nonlinear 
response can be important. In this talk, we study nonlinear response using 
transport theory quantum field theory, and explain the connection between 
these approaches. 

In a system that is out of equilibrium, the existence of gradients in ther- 
modynamic parameters give rise to thermodynamic forces which lead to de- 
viations from the equilibrium expectation value of the viscous shear stress: 

<5(^^.)=77<i)i/^. + r;(2)ijJ2^... (l) 

where is the four dimensional four- velocity field which satisfies 

— 1. The first coefficient is the usual shear viscosity. The sec- 
ond has has not been widely discussed in the literature - we will call it the 
quadratic shear viscous response. 

The Boltzmann equation can be used to calculate tranfiport properties 
for weak coupling A^"* theory with zero chemical potentials. We introduce 
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a phase space distribution function f{x,k) (the underhned momenta are on 
shell). The form of f{x,k) in local equilibrium is, 

/^°^ - ,,i.)u,l)k. _^ --=r^k; N,:=l + 2n,. (2) 

We expand / around f'^ using a gradient expansion in the local rest frame 
where u{x) — 0. We keep only linear terms that contain one power of H^,^ 
and quadratic terms that contain two powers of H^^,,. We write, / — + 
f'"^^ + f'"^^ + • • •• The viscous shear stress tensor is given by 

Using the gradient expansion of / to calculate {'Kij) and comparing with (|l|) 
we have, 



-15 7 (2^-(l + -)^ ^® 



We show that B(k) and C(fc) can be obtained from the first two equa- 
tions in the hierarchy of equations obtained from the Enskog expaiision of the 
Boltzmann equation. The first order equation can be cast into eH3, 

h,{k) = \ J^^^dT,2^,kdn[B,,{p^) + B,,{p^) - B,,{k) - B,,{p^) ] (6) 

where d- ~ (l + ni)(l + n2)"-3/(l + ^fe)- The second order Bolzmann equation 
leads toB, 

Nkhj{k)Bi„,{k) = 1 I drd„ {[iVlCy7m(pJ + A^2Cy7mbJ " A^feCy/,„(fc) 
^ Jl23 

-iVga.imbg)] + l[Ni2B,,{p^)Bi„,{p^) ~ Nk3B.,,{p^Bi^{k) + N:iiB,j{p^)Bi^{p^ 

+NklB,j{p^)Birn{k) + N32B,j{p^)Bi^{p^) + Nk2B,j{k)Blrn{p^)]} 

where we used Nij — Ni + Nj, Nij — Ni — Nj — 1, 2, 3, k). This equation 
can be solved self consistently for the quantity Cy7,„ (fc) using the result for 
B,,ik) from (I). 

Now we turn to response calculation □. We work with the den- 
sity matrix in the Heisenberg representation which satisfies §f = 
and write p = e"^+^/Tre-^+^ where A = J d^xF^To^ and B(t) = 
J d^x f_^dt'e''^^'-*^T^^ix,t')dt'F''{x,t') with = Pu^' and e to be taken 
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to zero at the end. Here A is the equihbrium part of the Hamiltonian and B 
is a perturbative contribution. We expand the density matrix in B and find 
the shear viscosity 

V^'^ = ^,:^M^n^Dnm\,o=o- (7) 
10 dqo q^t) 

This is the weh known Kubo formula i'S. The quadratic shear viscous result 
can be written as a retarded three-point correlator B: 

70 dqo dq'Q q^o " 

We have obtained a type of nonlinear Kubo formula that allows us to obtain 
the quadratic shear viscous response from a retarded three-point function 
using equilibrium quantum field theory. 

Next we obtain a perturbative expansion for the correlation func- 
tions of composite operators Dii{x,y) and GRi{x,y,z) which appear in 
(|) and §). We use the CTP formulation and work in the Keldysh rep- 
resentation. We define the vertices and Mijim by truncating exter- 
nal legs from the following connected vertices: Tfj = {TcTTij{x)(j){y)(j){z)) , 

^ijlni = {TcT^ij{x)TTlm{v)4l{z)(t){'w)) whcrC TTy (x) = di(j){x)djcj){x) - 

\5ij{dm<j){x)){dm<i}{xy). These definitions allow us to write the two- and three- 
point correlation functions as integrals of those vertices. Using the Kubo 
formulea above we obtainEI, 
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dKk^ PkUkil + rik 

dKk^Pknk{l + nk)Nk 



RerR2(K) 



ImEk 

ReMRi(K) 



ImSL 



(8) 

(9) 



Comparing with (^) and we see that the results are identical if we identify 



Bik) = 



RerR,2(k) 
ImSk 



Cik) = - 



ReMRi(k) 
ImSk 



(10) 



Hh 



) where m 



with the momentum K on the shifted mass shell: 
+ Resu- 
lt is well known that ladder diagrams give the dominate contributions 
to the vertex Tij. They contribute to the viscosity to the same order in 
perturbation theory as the bare one loop graph and thus need to be included 
in a resummation. The integral equation that one obtains from resumming 
ladder contributions to the three-point vertex has exactly the same form as the 
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equation obtained from the linearized Boltzmann emiation ^ with a shifted 
mass shell describing effective thermal exciiations WB. 

Following the pinch effect argumcntcrla, one can show that an infinite 
set of ladder graphs and some other contributions which we will call extended 
ladder graph contribute to the same order to vertex Mijim as the tree diagram. 
Therefore, for consistent calculation, we consider an integral equation which 
resums all of these diagrams, as shown in Fig. 1. 




+ 1/2 



(a) (b) 
Fig. 1: Integral equation for an extended-ladder resummation. 

This Integral equation can be cast intoi: 
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Note that once again we have obtained an integral equation that is decou- 
pled: it only involves Af/fi and Tji2. With V]^2 determined from the integral 
equation for the ladder resummation, can be solved to obtain Mm. Fi- 
nally, comparing (^) and (11) with (|^) and (|^) we see that calculating the 
quadratic shear viscous response using transport theory describing effective 
thermal excitations and keeping terms that are quadratic in the gradient of 
the four-velocity field in the expansion of the Boltzmann equation is equiv- 
alent to calculating the same response coefficient from quantum field theory 
at finite temperature using the next-to-linear response Kubo formula with a 
vertex given by a specific integral equation. This integral equation shows that 
the complete set of diagrams that need to be resummed includes the standard 
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ladder graphs, and an additional set of extended ladder graphs. Some of the 
diagrams that contribute to the viscosity are shown in Fig. 2. 








Fig. 2: Some of the ladder and extended ladder diagrams that con- 
tribute to quadratic shear viscous response. 



This result provides a diagrammatic interpretation of the Chapman- 
Enskog expansion of Boltzmann equation, up to quadratic order. 
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